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Abstract 

The first 2+ states in N=20 isotones including neutron-rich nuclei ^^Mg and '^^Ne are studied 
by the Hartree-Fock-Bogoliubov plus quasiparticle random phase approximation method based on 
the Green's function approach. The residual interaction between the quasiparticles is consistently 
derived from the hamiltonian density of Skyrme interactions with explicit velocity dependence. 
The B{E2,0f — > 2^) transition probabilities and the excitation energies of the first 2+ states are 
well described within a single framework. We conclude that pairing effects account largely for the 
anomalously large B(E2) value and the very low excitation energy in ^^Mg. 
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I. INTRODUCTION 



The pioneering observation in 1975 of the anomalous binding energy gain in very neutron 
rich Na isotopes revealed the breaking of the N=20 shell closure and the possibility of 
deformation The evidences of the breaking of N=20 shell closure in neutron-rich Mg 
and Ne isotopes are more clearly seen from the observations of E2 properties, the large 
B(E2]value in ^^Mg and the low excitation energies of the first 2+ states in ^'^Mg and 
30Ne 0. 

Several theoretical studies have been done to describe the anomalous binding energy and 
E2 properties in neutron-rich nuclei around N=20. Constrained Hartree-Fock (HF) calcu- 
lations of Na isotopes have been performed and ^^Na was suggested as deformed. Early 
studies made by Wildenthal and Chung showed that shell model calculations within the sd 
shell model space cannot explain the extra binding energies in this region. Subsequent shell 
model calculations |1, 0, IE 0] have demonstrated that the inclusion of the fp shell into the 
shell model active space is essential. The effects of the breaking of the N=20 shell closure are 
clearly shown in the description of the B(E2) values and the excitation energies in ^^Mg and 

30t 



ysnown in 

Ne 0, 11, 3|. The neutron 2p-2h configurations across the N=20 shell imply deformation 
of these nuclei. However, in the framework of the mean-field approximation with pairing cor- 
relations, like Skyrme Hartree-Fock-Bogoliubov (HFB) calculations, the calculated ground 
states in ^^Mg and ^''Ne turn out to be spherical (see, e.g., [11 Q). One possible way 
to describe the (dynamical) deformation is to include correlations beyond the mean field. 
Generator coordinate method (GCM)0, 3, 17 1 and antisymmetrized molecular dynamics 
calculations have been done in this direction. 

Nevertheless, the experimental evidence of deformation in ^^Mg is not well established. 
The energy ratios of the first 4+ and 2+ states, E{A'^)/E{2^), are 3.0 in ^^Mg and 3.2 in 
^"^Mg [31, and these values are undoubtedly close to the rigid rotor limit of 3.3. On the 
other hand, the ratio is 2.6 in ^^Mg, and this value is in between the rigid rotor limit and the 
harmonic vibration limit 2.0 [ilij. Moreover, the B(E2) value (in single-particle units) is 



15.0±2.5 in '^^Mg. This value is larger than in the other stable N=20 isotones but smaller 
than in other deformed Mg isotopes (21.0±5.8 in ^"^Mg, and 19.2±3.8 in ^^Mg). 

Generally speaking, the neutron 2p-2h configurations can originate not only from defor- 
mation effects but also from neutron pairing correlations. In the ^^Mg nucleus these two 
effects may coexist and make the large B(E2) value and the low excitation energy of the 2"*" 
state. In shell model studies it is not clear which effect is more essential to describe these 
anomalous properties. 

Recent angular-momentum projected GCM calculations with the Gogny force [H ^| 
were successful to reproduce the systematic trend of the B(E2) values and the excitation 
energies of the first 2+ and 4"*" states in Mg and Ne isotopes. However, the calculated 
excitation energies of the first 2+ states are somewhat higher than the experimental data 
(e.g., about 1.5 MeV in ^^Mg and about 2.1 MeV in ^°Ne). These discrepancies may be 
explained by the weakness of neutron pairing correlations around the spherical ground states 
in ^^Mg and ^°Ne in the corresponding HFB mean fields. 

The purpose of this paper is to emphasize how neutron pairing correlations play an es- 
sential role in the description of E2 properties in ^^Mg and ^°Ne. The existence of neutron 
pairing correlations means the breaking of the N=20 shell closure. As we will see, the ap- 
pearance of neutron pairing correlations is related to a special mechanism in loosely bound 
systems. We study the first 2+ states in N=20 isotones in the framework of self-consistent 
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quasiparticle random phase approximation (QRPA) with Skyrme interactions 21|. The 



QRPA equations are solved in coordinate space by using the Green's function method 2^ . 
Spherical symmetry is assumed for simplicity. The residual interaction between the quasi- 
particles is self-consistently derived from the hamiltonian density of Skyrme interaction that 
has an explicit velocity dependence. We will show that the B(E2) values and the excita- 
tion energies of the first 2"*" states in N=20 isotones, from the stable nucleus ^^Ar to the 
neutron-rich nuclei '^^Mg and ^''Ne are well described within a single framework and a fixed 
parameter set. The paper is organized as follows. In Sect. |n]we briefly describe the HFB 
plus QRPA calculations that we have done. In Sect. II I II we present the general results for 
the ground states of the N=20 isotones studied here. In Sect. |TV|we discuss the calculated 
and experimental E2 properties of these nuclei. Conclusions are drawn in Sect. 



II. HFB-QRPA CALCULATIONS 
A. Formulation 

We use the approach of self-consistent HFB-QRPA calculations with Skyrme interac- 
tions 21, 2^1. By self-consistent we mean that the HFB mean fields are determined self- 
consistently from an effective force and the residual interaction of the QRPA problem is 
derived from the same force. The QRPA problem is solved by the response function method 
in coordinate space. A detailed account of the method can be found in Ref.ji^]. Here, we 
just recall the main steps of the calculation. The QRPA Green's function G is solution of a 
Bethe-Salpeter equation, 

G = Go + GoVG . (1) 

The knowledge of G allows one to construct the response function of the system to a general 
external field, and the strength distribution of the transition operator corresponding to the 
chosen field is just proportional to the imaginary part of the response function. 
In Eq.(l) the unperturbed Green's function Gq is defined as 

° ^ ' ' ' ^ huj-{Ei + Ej)+ Z7] huj + {Ei + E^) +ir) ' ^ ' 
where the functions W{ra) are introduced as 

Wi,,{ra) = U,{ra)U,{ra) V,{ra)V,{ra) . (3) 
\-V,ira)V,ircx) -U,{ra)U,{ra) j 

Here, the U{r),V{r) are quasiparticle wave functions, the index a (a=l,2,3) stands for 
particle- hole (ph), particle-particle (pp) and hole- hole (hh) channels. The notation f{ra) = 
— 2(j/(r — cr) indicates time- reversal and = Wij — Wj^i. 

The residual interaction V between quasiparticles is derived from the Hamiltonian density 
< if > of Skyrme interaction by the so-called Landau procedure, 

14, {rar, r'a'r') = . . . • (4) 

dpi3 [r'a'T') dpa {rar) 
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The notation a means that whenever a is pp (hh) then a is hh (pp). The normal and 
abnormal densities are defined as 

< 0|^(m)V^(r(T)|0 > . (5) 
< Ol'ip^lra)^^ {ra)\0 > J 

The residual interaction V has an explicit momentum dependence, 

V(r, r') = f(Ku + "Av, '^u + ^v, "^u ± Vy, ± ^v]S{r - r'). (6) 



^ Pphircr) \ 






ppp{,ra) = 




H 


\phh{rcr) J 







The explicit form of the form factor F is shown in Ref. |21|. The operators with ^ (— act on 
the coordinate r (r'), and the operators with the index U (V) operate on the quasiparticle 
wave functions U{r) {V{r)) only. These momentum dependence are explicitly treated in our 
calculation. Because we calculate only natural parity (non spin-fiip) excitations, we drop the 
spin-spin part of the residual interaction. The Coulomb and spin-orbit residual interactions 
are also dropped. 

We have studied the influence of these momentum dependent terms and we have found 



that they can be important [21|. In a fully consistent calculation the spurious center-of- 
mass state should come out at zero energy. In practice, the full velocity dependence of the 
residual interaction (0) is replaced by a Landau-Migdal form to reduce the computational 
efforts (2^. Then, the self-consistency between the mean-field and the residual interaction 
is broken. To recover the self-consistency approximately, the residual interaction has to 
be renormalized by a factor adjusted so as to have the spurious stcitc cit zero energy. For 
example, the renormalization in Ref. 2^1 is about 20%. If the residual interaction in Eq.® is 



used, the renormalization factor is very close to 1 and the self-consistency is well recovered. 
For comparison, if one drops all momentum-dependent terms in Eq.® the renormalization 
factor would be about 0.6 . 

Another important aspect is related to the description of the low-lying states. The 
B(EA) transition probability is very sensitive to the treatment of the residual interaction. 
For example, the B{E2, O'l 2'^) in calculated with SkM* parameter with (without) 
these momentum dependent terms is 34.1 (20.9) e^fm^. The B(E2) value increases by 64 % 



21[. The experimental data is 28 ± 2 e fm [23[ which is close to the value calculated with 



the momentum dependent terms. Thus, the full residual interaction (jU)) is important for 
describing quantitatively the low-lying states and for comparison with experimental data. 



B. Inputs 

We apply the above formalism to study the first 2+ states in N=20 isotones, ^°Ne, ^^Mg, 
^^Si, ^^S and ^^Ar. The ground states are given by Skyrme-HFB calculations. The HFB 
equation is diagonalized on a Skyrme-HF basis calculated in coordinate space with a box 

Spherical symmetry is imposed on quasiparticle wave 

50 MeV, and the angular 



boundary condition 

functions. The quasiparticle cut-off energy is taken to be Ecut 
momentum cut-off is Imax = 7h in our HFB and QRPA calculations 



The Skyrme parameter sets SkM*[27| and SkP[28| are used for the HF mean-field, and 



the density-dependent, zero-range pairing interaction 



Vpair (^7 ^ 



^pair 




Sir- r') 



(7) 
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is adopted for the pairing field. The parameters a and pc are fixed as a = 1 and pc = 0.16 
fm~^. The strength Vpair is determined so as to reproduce the experimental neutron pairing 
gap in ^°Ne, An,exp.{^^^^) — 1-26 MeV. ^°Ne is the lightest mass even-even N=20 nucleus. 
The experimental pairing gaps are extracted by using the 3-point mass difference formula ,29|, 

A„(A^) = A'^\N - 1) = ^^[E(A^ - 2) + ^(A^) - 2E{N - 1)]. On the other hand, 
the average pairing gap in HFB calculations is defined as the integral of the pairing field, 
= J drpn{r) Anif) / J drpn{r) [s^. The pairing strength adopted for SkM* is Vpair = —418 
MeV fm~^, and for SkP, Vpair = —400 MeV fm~^. Fig Q shows the experimental and the 
calculated pairing gaps in 26,28,30]\j-g_ With these Skyrme parameters and pairing strengths, 
we get finite pairing gap in ^°Ne (vanishing of N=20 shell gap) and zero pairing gap in ^^Ne 
(appearance of N=16 shell gap) at the same time. 



III. GROUND STATE PROPERTIES 

FigO shows the neutron single-particle levels in N=20 isotones calculated in HF with 
the SkM* force. Results with the SkP force are qualitatively the same. The N=20 shell 
gaps change from 4.2 MeV in ^''Ca to 3.4 MeV in ^°Ne. The N=16 shell gaps change from 
2.4 MeV in ^°Ca to 4.0 MeV in 3°Ne. Within HF we can describe the vanishing of N=20 
magi city and the appearance of N=16 magic number at the same time. The change of N=20 
shell gaps looks moderate in comparison with the effective single-particle energies in shell 
model calculations 3. However, the definitions are different. In shell model calculations the 



single-particle energies are inputs of calculations and they are determined so as to reproduce 
the neutron separation energies and the one-particle spectra of ^^O and ^^Ca, and the change 
of the effective single-particle energies according to proton number are due to the change of 
the many-body correlations. On the other hand, the change of the single-particle energies 
in mean-field calculations reflects the self-consistent change of the mean-field potential. 

An important feature in FigO is the behavior of low-/ orbits, 2^3/2 and 2pi/2 in the 
fp shell. As the proton number decreases the single-particle energies of the high-/ orbit 
1/7/2 change almost linearly while the changes of 2p3/2 and 2pi/2 energies become very slow 
around zero energy. Moreover, the spin-orbit splitting of 2^3/2 and 2pi/2 states becomes 
smaller. As pointed out by Hamamoto et a/.jsH these effects can be understood by different 
/-dependences of the kinetic energy and the spin-orbit form factor as the single-particle 
energy comes close to zero. Because of these different /-dependences of the single-particle 
energies, the level density in the fp shell becomes higher with decreasing proton number, and 
the three orbits I/7/2, 2^3/2 and 2pi/2 become almost degenerate in ^"Ne. We can describe 
this behavior naturally by solving the HF and HFB equations in coordinate space but it is 
difficult to get this property by the methods based on the harmonic oscillator basis. 

FiglSl shows the HFB neutron and proton pairing gaps in N=20 isotones calculated with 
SkM* and SkP. The pairing strengths are adjusted so as to reproduce the experimental 
neutron pairing gap in ^''Ne. As the proton number increases, the neutron pairing gaps 
decrease monotonically and eventually, the neutron pairing gap becomes zero (for both SkM* 
and SkP) in ^^Ar as expected in stable N=20 nuclei. The interesting point is that the N=20 
shell gap itself changes very moderately but the calculated neutron pairing gap changes 
considerably from 1.26 MeV in ^°Ne to zero in ^^Ar. The mechanism can be understood by 
the increase of the level density in the fp shell when the proton number decreases, as noted 
above. Since the neutron pairing gap is adjusted in ^°Ne it remains close in ^^Mg for both 
interactions, but large differences are seen in ^^Si and '^^S. On the other hand, the calculated 
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proton pairing gaps do not depend on the Skyrme force. 

Fig|3 shows the average number of neutrons Nfp in the fp shell in N=20 isotones, cal- 
culated in HFB with SkM* and SkP. According to the change of the neutron pairing gaps, 
Nfp decreases monotonically from ~ 0.8 in ^°Ne to ~ 0.5 in ^^Mg. These values are very 
different from the prediction of the " island of inversion" , Nfp = 2 |9| and Monte Carlo shell 
model, Nfp > 2 in ^ONe and ^^Mg [ij. 

IV. B(E2) VALUES AND EXCITATION ENERGIES 

We have calculated the first 2+ states in N=20 isotones in HFB-QRPA with Skyrme 
interactions, assuming spherical symmetry. At the mean-field level the ground states in 
N=20 isotones including ^^Mg and ^°Ne have been found to be spherical (see, e.g.,, jl3l.E3|). 
Our aim is to investigate whether these states can be described as vibrational states built 
on the spherical ground states. 

In Figl^lour results of QRPA calculations with SkM* are compared with the predictions 
of the Monte Carlo shell model (MCSM) and the available experimental data 0, 0, 
0, 113, 13 ■ Our QRPA results are in good agreement with the experimental data and 
they are consistent with the MCSM prediction of the B(E2) value in ^°Ne which is not yet 
measured experimentally. The QRPA calculations have been done with the SkM* parameter 
set and the fixed pairing strength Vq = —418 MeV fm~^ the choice of which is explained in 
subsect. 2.2 . The general properties of the first 2"^ states in N=20 isotones, namely very 
large B(E2) values and very low excitation energies in ^^Mg and ^°Ne are well reproduced 
within a single framework. 

To check the interaction dependence we have carried out QRPA calculations with the 
SkP interaction. FigE] shows the B(E2) values and excitation energies of the first 2"*" states 
with SkM* and SkP. Concerning the B(E2) values we get similar results, especially very 
large B(E2) values in ^^Mg and ^°Ne. On the other hand, large differences are seen in 
the excitation energies in ^^Si and ^^S. This can be understood by the difference in the 
neutron pairing correlations shown in Fig|H| In '^°Ne, '^^Mg and ^^Ar the neutron pairing 
gaps calculated with SkM* and SkP are almost the same while they are somewhat different 
in ^^Si and ^^S. Because the neutron pairing gaps are larger in SkP than in SkM*, the 
excitation energies become lower with SkP than with SkM*. 

We now explain how the neutron pairing correlations are important to make the B(E2) 
values larger and the excitation energies lower. To see which two-quasiparticle configurations 
contribute to make the low-lying 2+ states, we show the unperturbed isoscalar quadrupole 
strength functions in N=20 isotones calculated with SkM* in FigjTI The peaks indicated by 
straight (dotted) arrows correspond to proton (neutron) two-quasiparticle configurations. 
All these neutron two-quasiparticle configurations appear because of the neutron pairing 
correlations. Many peaks of the neutron configurations are seen in 3°Ne, ^^Mg. On the 
other hand, the neutron configurations are negligible in ^^Si and they completely disappear 
in ^^S. The B(E2) values are primarily made of the proton configurations in the sd-shell 
but the neutron configurations assist to make the B(E2) values larger by coherence between 
protons and neutrons. Actually, if the neutron pairing correlations are neglected the B(E2) 
values become very small and the excitation energies are sizably higher in '^^Mg and ^''Ne, 
as shown in Fig|Sl Under these considerations, we can conclude that the very large B(E2) 
values and the very low excitation energies in 32 Mg and ^"Ne appear thanks to the presence 
of the neutron pairing correlations. We have seen in sect. Illll that. around the drip line the 
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origin of neutron pairing correlations lies in the different behavior of the single-particle levels 
with different orbital angular momentum / as the levels approach the separation threshold. 

Generally speaking, neutron 2p-2h configurations across N=20 can originate from defor- 
mation effects or pairing effects. Both effects can a priori contribute in the nucleus '^^Mg. 
We have shown that a spherical QRPA description, i.e., putting emphasis on the pairing 
aspects and neglecting the possible deformation effects, can give very satisfactory results. In 
the previous studies based on shell model calculations 3, U) ^| the importance of neutron 
2p-2h configurations for describing the B(E2) values in ^^Mg, ^°Ne was shown, but the re- 
spective roles of pairing and deformation were not clear. The angular momentum projected 
GCM calculationsHi, |l3 are successful in reproducing the general trend of the B(E2) val- 
ues and the excitation energies in Mg and Ne isotopes. However, their predicted excitation 
energies are somewhat higher than experiment. This may be due to the weakness of neutron 
pairing correlations in these calculations. 

Before closing this section we would like to make a brief comment on the stability of 
QRPA solutions. As it is well known, when a QRPA eigenvalue is approaching zero the 
solution suffers instability and the transition probability diverges. Since our QRPA solutions 
for ^^Mg and ^'^Ne have very low energies, we have to check whether the calculated B(E2) 
values are really meaningful or just spurious results. FigEl shows the pairing strength Vpair 
dependence of the excitation energy and the B(E2) value in ^^Mg. If \Vpair\ increases, the 
excitation energy decreases and the B(E2) value increases. This behavior is the result of 
competition between two effects. First, the pairing gap and also the quasiparticle energies 
increase with increasing |V^air|- Therefore, the two-quasiparticle energies and the QRPA 
excitation energies should increase and the B(E2) values should decrease. Second, if \Vpair\ 
increases, many two-quasiparticle configurations can contribute to make the 2"*" state and 
the collectivity increases. In this case the QRPA excitation energies would decrease and the 
B(E2) values would increase. In the ^^Mg case the second effect wins (cf. FigEl). As \Vpair\ 
is increasing, the excitation energy becomes lower and the B(E2) value increases linearly 
up to \Vpair\ ^ 422 MeV fm'^. Above \Vpa^r\ ^ 426 MeV fm'^ the B(E2) value starts to 
diverge. Because our adopted pairing strength is Vpair = —418 MeV fm~^, we confirm that 
our QRPA solution is meaningful. 



V. CONCLUSION 

We have studied the first 2+ states in N=20 isotones by the HFB-QRPA model with 
Skyrme interactions. The residual interaction is consistently derived from the Skyrme Hamil- 
tonian, keeping all its explicit momentum dependence. 

Because of the different behaviors of the neutron I/7/2, 2^3/2 and 2pi/2 levels when the 
single-particle energies are approaching zero, the neutron pairing gaps have finite values. 
This mechanism breaks the N=20 magicity in ^^Mg and ^°Ne. 

Within our consistent QRPA calculation with spherical symmetry the B{E2,0^ 2^) 
values and the excitation energies of the first 2"*" states in N=20 isotones including '^^Mg 
and ^"Ne are well described. The existing experimental data are reproduced quantitatively. 
The B(E2) value in ^°Ne has not been measured yet but the QRPA value is consistent 
with the prediction of the MCSM. The important role of the neutron pairing correlation is 
emphasized. If the neutron pairing is dropped, we cannot get the correct B(E2) value and 
excitation energy in ^^Mg and ^°Ne. In the real ^^Mg nucleus, both neutron pairing and 
deformation effects may coexist and help to make the large B(E2) value, but our calculation 
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shows that neutron pairing correlations are essential. 
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FIG. 1: HFB neutron pairing gaps in 26, 28,30 j^^g calculated with SkM* and SkP. The pairing 
strengths Vpair are fixed so as to reproduce the experimental neutron gap in ^°Ne. The exper- 
imental pairing gaps are extracted by using the 3-point mass difference formula |29| . 




FIG. 2: HF neutron single-particle levels in N=20 isotones calculated with SkM*. Solid lines 
correspond to bound and resonance-like states, dashed lines to positive energy discretized states. 
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FIG. 3: The neutron and proton pairing gaps in N=20 isotones calculated in HFB with SkM* and 
SkP. 




FIG. 4: Average number of neutrons in the fp shell in N=20 isotones calculated in HFB with 
SkM* and SkP. 
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FIG. 5: The B{E2,0f — > 2f) transition probabilities and excitation energies of the first 2"*" states 
in N=20 isotones calculated in QRPA with SkM*. For comparison the predictions of MCSmQ 
and the available experimental data 0, IE Eli 2^ are shown. 
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FIG. 6: The B{E2,0f — > 2^) transition probabilities and excitation energies of the first 2'^ states 
calculated in QRPA with SkM* and SkP interactions. 
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FIG. 7: Tlie unperturbed isoscalar quadrupole strength functions in N=20 isotones calculated 
with SkM*. The peaks indicated by straight (dotted) arrows correspond to proton (neutron) 
two-quasiparticle configurations. 
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FIG. 8: The B{E2,0^ 2^) values and the excitation energies of the first 2+ states in N=20 
isotones calculated with/without neutron pairing correlations. Proton pairing is included in both 
cases. 
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FIG. 9: The B(E2) value and excitation energy of the first 2"^ state in ^^Mg calculated in QRPA 
with SkM*, as a function of the pairing strength Vpair- 
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